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O . Abstract 



We present an interesting connection between Einstein-Friedmann equations for the models 
of universe filled with scalar field and the special form of Abel equation of the first kind. This 
connection works in both ways: first, we show how, knowing the general solution of the Abel 
equation (corresponding to the given scalar field potential) one can obtain the general solution of 
the Friedman Equation (and use the former for studying such problems as existence of inflation 
with exit for particular models). On the other hand, one can invert the procedure and construct 
the Backlund auto-transformations for the Abel equation. 
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1. Introduction 



(N 
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Starting out from the classical work of Abel [1] , the Abel equation has seen a lot of extensive studies 
and applications [2], [3] (see also [4], [5]). The reason for later lies in the fact that Abel equations 
rather frequently appear in the process of reduction of order for many second (and higher) order 
families [6], [7], and hence are often found in the modeling of real problems in varied areas. 

In this work we present but another new field of application of the Abel equation: as we shall 
see, the general solution of the cosmological Einstein-Friedmann equations for the universe filled 
with scalar field for a given potential can be expressed via the general solution of the Abel equation 
of the first kind. 

Let us start by writing the Friedmann equations describing homogeneous isotropic universe (in 



Friedmann-Lemaitre- Robertson- Walker metric) : 

• dV 

4> + m4> + -j- = o 

i # k w 

H 2 = -4> 2 + v - - 
I o z 

where a is a scale factor, cj) - scalar field, V is the scalar field's self-action potential, H - Hubble 
constant H = a/a, k = 0, ±1 and we have incorporated the rescaling = c = 1. 

The general study of the dynamics of the universe for a given self-action potential is an excep- 
tionally difficult mathematical problem. Nevertheless, this paper is dedicated to the task of finding 
out a regular procedure of construction of the system ([T])'s general solution for a given potential. 
By definition, the general solution has to have two arbitrary integration constants: = cf>(t; to, cjE 



1 Constant t is connected to the translation invariance of {TJ 
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Restricting ourselves to the case of flat space-time k = 0, it would be enough to define the general 
solution for 0, since for such model it is possible to exclude the Hubble constant from the first 
equation of ([!]) and rewrite it as 



Scale factor a(t;to,C) then would be defined by simple integration of the second equation of (pQ), 
the process that will give birth to the second (multiplicative) integration constant, allowing for the 
scale factor rescaling. 

Remark 1. For any twice differentiable scale factor function a = a(t) and for any k it is possible to 
find out the corresponding scalar field that will provide the necessary dynamics. This can be done 
via the following formulas: 



12 f , y/a 2 — a'd + k 
(t) = 0o ± \ 77 / dt 



V(t) 



3 . 

9 , ( 2 ) 
ad + 2d + 2k 

3^2 • 



Equations ((2J) describe the function V = V((f>) in the parametric form defined for any given function 
ait). Of course, this procedure has nothing to do with the searches for the general solution of 
JT]). First of all, the described procedure might at best help to introduce the first of the constants 
(4>o, obtained via the integration of the first equation of (J2j) ) , but not the second oneH It's second 
(and worst from the physical point of view) fault lies in the fact that the potential has to be 
derived from the known scale factor, and not the other way around - in this sense, it can be called 
the "inverse" problem. As easy as the "inverse" problem goes, the "direct" problem is way more 
difficult - it would be enough to say that in most cases it is being solved only by the procedures of 
the numeric calculations. However, as we will show below, it is possible to reduce it to one of the 
most well-developed problem: to finding the solution of the particular Abel equation of the first 
kind. 

2. The functional of full energy 

In the rest of the article we'll restrict ourselves to the case k = as the one best fit to describe 
the observed universe. Following the ideas of the [8], let us introduce the functional of full energy 
(hamiltonian) W: 

W=^ 2 + V{4>). (3) 
It is easy to see that, using ([3]) and assuming k = system (JTJ) can be rewritten as: 

# (4) 
H = ±VW 



2 It is by all means possible to choose the scale factor depending on the arbitrary amount of undefined constants 
{Ci, Cjv} so that they will appear in the solution cj> — <f>(t, Ci, Cjv). But then, as follows from ([2]), the potential 
V will depend on those constants as well. 
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Moreover, knowledge of the function W = W(4>) provides an elementary way of finding all other 
quantities ((f) = (f>(t), a = a(t), V = V ((/))). If W ^ (cf. Remark 2), <f>(t) can be derived from the 
ordinary differential equation 

.1 w'{<t>) 



dt T 3 v /Wl>) 



and potential V: 



v(<j>) = w{4>) 



1 {W'i^f 



(5) 



(6) 



18 W((/>) ' 

As for the scale factor a(t) = exp (J H(t)dt), it can be obtained by substituting (f>(t) from ([5]) into 
(j3J) and consequent integration. 



Example 1. Choose W(4>) = A</> 4 /4, A > 0. According to 
to the potential with the spontaneously broken symmetrjH 



V(cf>) 

whereas the dynamical variables would be: 



\<P 4 
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), this hamiltonian corresponds 



(7) 




4v / A(t-t ) 



(t - 1 ) 



(8) 



VA^ ( , 4>/A(t-t j 



Here the variables with zero subscript correspond to current (t = to) values of scalar field and scale 
factor. 

Integration of © gives rise to the constant to, which appears because of the translational 
invariance of (pQ) and can be assimilated by the translation t — > t + const. In order for the solution 
(p = 4>(t;to,C) to be general it has to contain the second independent constant C. This constant 
can be obtained through the following calculations: consider |6} as a differential equation w.r.t. 
variable W((p) with the given V((f>). General solution of this equation contains the thought after 
integration constant C . Thus, we can now formulate the following 

Proposition 1 If for a given V(4>) the general solution of equation isW = W(4>, C), the general 
solution of the Friedmann equation (OJ) (p(t;to,C) exactly corresponds to the general solution of (El). 

Hence, the problem is reduced to the task of finding the general solution of © for a given V((f>). 
Remark 2. The case W = shall be considered separately. It is easy to verify that for any given 
V < one gets H = (a(t) = ao = const - stationary universe) and the general solution of the ([1]) 
contains the single integration constant: 



d((> 



±(t-t ); V^O 



V = 



0) 



'Note, that the value of constant fi in (0 isn't arbitrary: /j 2 = 4A/9. 
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So this model is meaningful for the non-positive potential V(4>) only [9]. It is also interesting to 
note that the parameter of equation of state w = p/p = oo where p is the density and p is the 
pressure of scalar field. 

In what follows we'll restrict ourselves to the case W 7^ 0. 

3. Main Theorem 

The main result of this paper lies in the following theorem: 

Theorem 1 Let x = 3\/2^>, x = l n |V|, k = ±1. For a given V(cf>) the corresponding hamiltonian 
W = W(x,C) is defined as (cf. Remark 2): 



W(x,C) = V(x) 



( 



y + vV-1 +1 



1- (y + Vv 2 -^ 2 



(10) 



where y = y{x, C) 7^ ±1 is a general solution of Abel equation of 1st kind: 

y' = -\{y 2 -i)^-x'y)- (n) 

Moreover, the special case V = occurs if and only if y = ±1 and the hamiltonian W has the form: 

W = Ce KX (12) 
The proof of the theorem can be performed by the direct calculations. 



Remark 3. (|10p defines a family of solutions of (pQ), parameterized by the constant C . Substituting 
W(x, C) in ([5]) after the integration one will obtain <p = <p(t; C, to), where to is the second integration 
constant, connected to the invariance of the scalar field relative to translations t — * t + const. In 
other words, the suggested algorithm indeed allows to find the general solution of ([1]), which is a 
main result of the paper. 

Remark 4. For the special case y = ±1 (i.e. V = 0) the general solution of the ([1]) has the 
form 

<t>{t) = 0o ± -y- log (t - t ) , 

the scale factor is 

a{t) = a (t- t ) 1/3 
and the case corresponds to the "stiff" equation of state with w = 1. 

Remark 5. If x' const (cf. Ch. 4, example D) the eq. (|11|) has two fixed points y = ±1. 
It means that for any initial yo = y(xo,C) 7^ ±1 the solution y(x,C) 7^ ±1 for all values of x 
excluding probably x = ±00. 

The general form of Abel equation of first kind is: 

3 

y' = E fn v"* ( 13 ) 

r?=0 
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where, in our case: 



(14) 



As known (cf., for example [5], [6]), if /i is continuous, /2 and are continuously differentiable 
and /3 / then one can represent the equation (jlip in normal form: 



rf = rf + J(x), 



(15) 



where 



y = w(x)t/(£) + 



u{x) 



: CXp 



1 fV'(x)oj 2 (x) 
V{x) 



dx. 



y/Vjx)~" r \ 6 7 V'(z) 
2k (W" - V) V 2 



dz 



(16) 



J 



(W'u) 



Example 2. For the popular cosmological model with quadratic potential 3 

, , m 2 d> 2 

one can use the substitution (fl~6l) for the <\> > (x > 0). In this case 



(17) 



n 6 . 

Loix) = e 

mx 



-x 2 /24 



2m 2 



J(x) 



Ei 1, 



12 



K1JI' J 

23328' 
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-x I x 



2 -18) e* 2 / 8 , 



-x 2 /12 



(18) 



Remark 6. In the case of polynomial potentials 



V 



n 



Xx r 



18«/ 2 



n 



with positive coupling A > 0, the Abel equation of 1st kind can be transformed into the particular 
case of the Abel equation of 2st kind. In this case one gets x' = n/x. Assuming |yo| > 1 and 
sign(yo) = —«, lets consider y as the independent variables and x = x(y) as a solution under the 
question. Introducing new function P = P(y) = nx(y) — ny and substituting it into the (fTT|) one 
get 

PP' = F 1 (y)P + F (y), (19) 



with 



Lets define 



Fi(y) 



-n 



F (y) 



2ny 



P = u{y)+F(y), F(y) = F F 1 {z)dz = -ny + log 



y-i 



4 In a quantum field theory this model describes noninteracting massive scalar particles. We shall carefully consider 
this model in the Sec. 5 
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Then the equation (|T9j) will be reduced to 

{u + F)u' = F . (20) 

Finally, if we'll assume that \yo\ > 1 then Fq ^ and one can introduce a new independent variable 

rv 

£ = / ^oO^tfa = -nlog (y 2 - l) , u(y) = 77(C)- 
Jo 

with whom the equation (|2U|) will has the normal form 

(n + F)rf = 1. 

Another form of the Abel equation of 1st kind (jlip can be obtained if we know at least one of 
it's exact solutions. For there are two such solutions y = ±1 = K. Calculating the function 

E{x) = exp J (3f 3 K 2 + 2f 2 K + fx) dx = V(x)e 

and using the substitution 



kKx 



y = K + 



E{x) 
z(x) ' 



(21) 



one get 



where 



For the case 



$2 = and the equation (|21|) can is exactly solvable. This example will be thoroughly examined in 
the next section. 



z' + — + $ 2 = 0, 

z 

= - (V 2 )' e- 2KKx , $ 2 = h~ KKx (3KV - kV) 
V = V e KKx/3 



Remark 7. There have been many works recently that treat either modified gravity or f(R)- 
gravity |10| . [llj . |12j . [13] . |14j . All these models suggest an alternatives for the origin of dark 
energy. It may be naturally expected that gravitational action contains some extra terms which 
became relevant recently due to the significant decrease of the universe curvature. 
The f(R) models can in general be describes via the action 

S =J?f d A x^-g[R + f{R)]i 

where f(R) is the proposed additional term. After the introduction of the Friedmann-Lemaitre- 
Robertson- Walker metric 

3 

ds 2 = -dt 2 + a 2 {t)Y J (dx i f 1 (22) 
i=i 

one gets new equations, naturally differing from the generic Friedmann equations (pQ). Clearly, the 
general solutions problem in this new framework gets even more difficult then before. 
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However, one can still reduce any given f(R) model to a special usual Einstein- Friedmann 
universe filled with scalar field with non minimal coupling. In fact, letting a(t) be the particular 
solution of some /(i?)-model in metric (|22p and substituting it into the ([2]) one will obtain the 
function V = V(4>) in the parametric form defined for the given scalar factor. Therefore, we end up 
with the following 

Proposition 2 For any solution a(t) of any f{R)-model in metric fljj there exist an Einstein- 
Friedmann model with scalar field having exactly the same scale factor. 

In other words, one can reduce /(-R)-models to usual Friedmann equation (JTJ) in order to find 
their general solutions. 

4. Some Examples 

Let's consider a couple of examples of potentials that are frequently used in cosmological studies. 



n 

Here we have the following set of functions f^: 

h = ~fi = h = ~fo = 



b) 



4 2 

This potential describes the popular field model of scalar particles with the coupling constant A > 
and the mass m. 

_ 2(Ax 2 + 18m 2 ) 

x(Xx 2 + 36m 2 ) 

c) 

V = A = const. 
h = ~fi = o. 

if 4> = const this case corresponds to the model with the cosmological constant. 



d) 



V = V e 6V * a *, a = const. 
h = ~fi = ol. 



Let us now consider the cases c) and d) in details. 
If V = A = const, equation (fTTj) has the solution: 



6 k(x-x ) t ^ 

y( x ' x o) = ^(x-xo)-! ' (23) 

and the functional PF has the form: 

W(x, x ) = A cosh 2 {^(x - x )) • (24) 
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After substitution into © we obtain the following differential equation: 

3k 



±kv'2A sinli ( ^=(0 - (f>o) 



(25) 



The general solution of (I25D (j)(x;xo,to) will be parameterized by two arbitrary constants xq, to and 
will have a form: 

\/2 

<j) = 0o + ~ — arccosh ( cot 



3 k 



The scale factor 



a(t) = clq sinh ( 3a/A 1 1 — to | 



3\/A|t - t 1 

±1/3 



(26) 



(27) 



Note, that if t — > oo then <j> — > (see (|26p ) so the density p — ► A and the pressure p — > —A. In 
the case of positive power in (|27p we obtain the well-known de-Sitter (dS) solution as t — > oo: 

<•(') - 

while the negative power case contribute the solution with the so called Big Rip singularity at t = to 
which means that one deals with the phantom cosmology here [15], [16], [17j . [18] . 
For the case d) equation (lllj) takes the form 



y' = a(y - i)(y + i)(y - s), 



(28) 



where s = n/(2a). This equation has three fixed points: y = ±1 and y = s. For simplicity one 
choose s > 1. Therefore if the initial value 1 < yo < s then this will be the case for y{x,C) for any 
x. The solution of the (12811 has the form 



y + 1 V (y ~ s ) 



s y-l) y 2 -l 
where C > is an integration constant. One can see that 



(j^2a(s 2 — l)x 



(29) 



for a < and 



y-> s 



y^l 



at 
at 



at 
at 



a; — > +oo 



-oo 



X 
X 



-oo 
-oo 



for a > 0. 



5. The "y model: inflation and slow-rolling approximation 

The Abel representation (jlip of Friedmann equation ([1]) can be extremely useful even in the cases 
when one cannot find it's exact solution. To demonstrate this let us consider the popular cosmo- 
logical model with the quadratic potential ()1T[) . It is known that Friedmann equation with this 
potential is non-integrable. Same goes for the Abel equation. However, by studying this model it 
is still possible to draw some interesting conclusions, namely: that the equation (pQ) may result in 
inflation with natural exit from it. 



S 



The standard approach to this task (cf., for example, [19]) would be to use the slow-rolling 
approximation, i.e. assume that 

K = f «m (30) 

If this is true, the equation of state is p ~ —p and one gets inflation. For polynomial potential 
(akin to the one in question) condition (|30p may be valid if (f) 3> 1. It is then assumed that during 
the inflation process the kinetic term K increases until the slow-rolling approximation (|30p is no 
longer applicable, which heralds the (spontaneous) natural exit from the inflation. This reasoning, 
being rather simple as an idea proved to be quite a challenge when it came to the rigorous proofs in 
particular cases. In the last 14 years this difficult problem has drawn a lot of interest and attention. 
The investigations conducted in these field (see [20] , [21] , [22] , [M] and list of references) have led 
to the following conclusions: 

1. The idea of slow rolling is true. Inflation does indeed occur under an extremely broad range of 
self-acting potential, and there is hence no need to fix a certain form of the potential to obtain an 
inflationary universe. 

2. The exit from inflation, on the contrary, turned out to be a problem. For many model potentials 
the universe never stops inflating. In general, the exit has to be achieved only by the means of 
fine-tuning or, in other words, by the parameter fitting. 

In this Section, by using the Abel equation we show that the popular non-integrable model 
with quadratic potential (|17p is free from the later problem, which lies in total accordance with the 
conclusions achieved so far. 

Let us start by introducing the quantity 9 2 (y) (cf. (fl~0|) ): 

I (y + ^Jy 2 -!) 2 + 1 ' 



(31) 



\l- (y + vV-i^ 
For y > 1 this is the monotonously decreasing function of y, such that: 

lim6> 2 (y) = +oo, lim 9 2 {y) = 1. 

V— >1 V— »oo 

The plot of 9 2 (y) is represented on the Fig.l. The slow-rolling approximation (|30p can be rewritten 
in the form 

^ = 9 2 {y) - 1 « 1. (32) 

For the y — > oo we have 9 2 {y) ~ 1 + 1/y 2 therefore the approximation (|32p will be valid if y 1. 

The inflation take place whenever a{t)/a{t) > (which is identical to condition p + 3p < 0). 
This will be the case if 9 2 (y) < 3/2 or y > \/3- The pressure will be negative if 9 2 {y) < 2 or y > 
All these results are presented in Tabl. 1. 

Now consider the Abel equation (TlTj) with the quadratic potential (fT7|) . i.e. with % = 2/x. For 
simplicity let us choose k = +1. The equation (|lip takes the form: 

<- -\ t" 2 - 1 ' ('-?)■ (33) 

Remark 8. Since y = ±1 are fixed points of the equation (|33|) the range of an arbitrary solution 
of ([33]) y(x,C)E will be splitted into the three distinct regions: h = {y \y > +1}; h = {y \ — 1 < 

5 In this chapter for the sake of simplicity we'll refer to the solutions of (|33p as just y(x), omitting the constant C. 
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y 



Figure 1: The plot of 6 2 (y). 





Slow-rolling < |V|) 


Inflation (p + 3p < 0) 


Negative pressure 


I: 1 <C y* < y < oo 


yes 


yes 


yes 


II: \/3 < y < y* 


no 


yes 


yes 


III: y/2 < y < y/3 


no 


no 


yes 


IV: y < V2 


no 


no 


no 



Table 1: 



y < +1}; and I3 = {y\y < —1}. Curve y(x) can cross the lines y = ±1 in two possible cases: (i) at 
x = 0, which is a singular point of (|33|) and (ii) if y(x) is singular in at least one point x = x s . In 
the framework of the discussed applicability of the Abel equations to the cosmological models with 
real scalar field, the most interesting regions would be I\ and I3. In this work we'll concentrate on 
studying I\ . 

Before we actually start tackling the equation (|33h . it would be worthwhile to spend some time 
on analyzing the easier but still somewhat similar Riccatti equation : 



(34) 



Similarly to (j33[) . equation (|34p also has a fixed point at y = 1 as well as the singular point at x = 0. 
However, contrary to ([33]) . this equation can be solved explicitly; it's general solution has a form: 



_ x + 2 + Ce*/ 2 



(35) 



where C is an integration constant. (|35p is obviously regular whenever C > 0, and C < implies 
existence of singularity at a certain point x s (see Fig. 2). Existence and location of singularity x s 
is uniquely determined by the initial conditions {a?o,yo} as follows: 



x - 2 log 



~ 2y - x - 2 ' 
. 2(yo-l) . 



(36) 
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Figure 2: The plot of solution (|35p at c = 1 (left) and c = — 1 (right). Horizontal line corresponds 
to the fixed point y = 1. Note, that y(:c) intersects it at x = 0. 



If x s , defined via (|36p is real, soultion (|35p will have a singularity there. In particular, choice of 
initial conditions lying in I\ and such that yo 3> xo will result in singularity (Fig. 2, right) at 
x s ~ xq. We would also like to point out that this choice in general leads to the violation of Lipshitz 
condition that in our case has the form 

\2(y 2 + yi) - x - 2\ < 2x L, 

where L is a positive constant. 

Note that in order to obtain the aforementioned estimate one can actually avoid the usage of 
solution ([35]) . but instead start from assuming that y — > +oo at s -> x s > and integrate the 
reduced version of ([Ml) y' = y 2 /x s to get 

y(x ~ x s ) = - Xs _ - . (37) 

(|37p can be obtained from (|35p with C = —2e~ Xs ^ 2 by the Taylor expansion around x = x s . Finally, 
the usage of the initial conditions {xo,yo} together with (I37p results in a crude estimate for x s (cf. 
(|36p . that gets more and more accurate for yo ^> 1: 

x , (38) 



yo - 1 



in a good accordance with the results obtained by studying the exact solution (|35p . 

As we have mentioned already, equation (|33p is not the one allowing for explicit integration, 
hence the study for existence and location of possible singularities has to be performed using the 
approximate approach that we have tested on the Riccatti equation (|34p . The analog of (|37p for 
(1331) has the form 



i.e. unlike the discussed singularities of solutions of (|35p there might exist at most one left-sided 
discontinuity point of second order, while for x > x s the solution is either finite or complex (cf., for 
example, solution of a simple ODE y' = y 3 ). By analogy with (|38p we'll obtain for (|33[) 

_ 2x y 2 
Xs ~ 2y 2 - 1 X °' 
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when yo ^> 1. Note that this condition together with y$ 3> 1 are also violating the Lipshitz 
condition, therefore it is expected for the singular point to appear near the critical point x s ~ xq. 
Contrary, the choice xq 1, a yo = 1 + o(l) is deemed to render the nonsingular solution (at least 
on the right half axes), by analogy with Fig. 2, left. For what follows we'll assume that the initial 
conditions are chosen in such a fashion that ensures continuity of solution of ([33]) at x > 0. 

Proposition 3 Solution of \33\) . continuous at x ^ is infinitely differentiable there. 

The proof of the proposition. Continuity of y(x) is postulated on the chosen class of initial 
conditions (cf. [23]). Differentiability of y(x) ar i / follows directly from (|33p . After n — 1 
differentiation of ()33f) w.r.t. x one gets the following: 



where k n = 2 n + l, and coefficients a^(x) are continuous and differentiable at x / 0, which concludes 
the proof. 

Theorem 2 Solution of (E3|) that belongs to I\ has one and only one maximum on the right half 
plane x > 0. 

The proof of the theorem. 

First, let us show that existence of extremum of solution y(x) on x > implies the uniqueness 
of maximum. 

Indeed, assume that x = x m > is an extremum of y(x), i.e. y'(x m ) = 0. According to ([33]) this 
means that either y = ±1 or y m = y(x m ) = x m /2. Since the second derivative at x m is equal to 



and x m > implies y m > it follows that inside of the I\ region y"(x m ) < 0, hence x m can only 
be the maximum. Uniqueness of maximum then follows straight from the analyticity of y(x) (cf. 
Proposition 3). 

Now consider the line y\{x) = x/2. According to all the above-said, if y{x) and yi(x) intersects 
at any point x m > 0, at this y{x) suffers a maximum and, moreover, this maximum is unique. Let 
us now prove that existence of such intersection point. Indeed, let us assume by contradiction, that 
y(x) and yi(x) does not intersect at x > 0. This means either one of two possibilities: 

A. y(x) < y^x) at x > or 

B. y(x) > y^x) at x > 0. 

Case A shall be excluded due to the fact that whenever < x < 2 we have y(x) > 1 (region 
but yi{x) < 1. Therefore, we are left with the case B, corresponding to the unbounded increase of 
function y(x). 

It is easy to see that function y(x) has no oblique asymptote at x > 0. Indeed, existence of 
such asymptote implies the existence of limit of left hand side of (|33p . corresponding to the angular 

6 Note that when x — > the solutions of (|33|) has a behaves as d n y/dx n — > - a very interesting property, differing 
it from solution (|35[) of Riccatti equation (|34[) . 
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coefficient of the asymptote, which contradicts the x 2 divergence of the right hand side of (|33p . 
Thus, y(x) has to increase fast enough, so that for large values of x the contributions of unities 
in both parentheses in (|33p will be negligibly small, which, in turn, leads to a simple differential 
equation y' = y 3 /x. General solution of this equation contains an unavoidable singularity at a 
certain point x s (C) > 0, whose location depends on the arbitrary integration constant, and for 
x > x s (C) the equation has no real solutions. □ 

This theorem allows for examinations of the asymptotic behavior of y(x) at x — > +00 and at 
x — ► +0 provided that y(x) belongs to the region I\. 

Proposition 4 For y{x) from I\ for x > the following limits exists: 

lim y(x) = 1, (40) 

x — > + oo 

lim y(x) = 1, (41) 

The proof of the proposition. The (lIOl) is a consequence of the following three facts: (a) the 
function y{x) has a maximum at x > 0; (b) the function y(x) belongs to I\ and (c) y = +1 is a 
fixed point for equation (f33l) . 

The limit (|4ip is somewhat more interesting. First, note, that despite the fact that y = 1 is a 
fixed point, since x = is a singular point of (|33p (by analogy with the Riccatti Equation, cf. Fig. 
2, left), (|4*T|) doesn't lead to a contradiction. In order to establish (|4*T|) it should be kept in mind 
that, according to the Theorem 2, y{x) is a bounded function, hence for sufficiently small values of 
x it is possible to neglect the unity at the right parenthesis of (|33p in order to get the ODE: 

, (y 2 - i)v 
y = 



2.r 



whose general solution is 



y(x) 



cx 2 



VT 

The (|4ip follows immediately. 

All these results can also be justified by the numerical integration; the results of one can be seen 
on Fig. 3. 

The initial value of field x 3> 1 so y ~ 1 and we have the "stiff" equation of state. During the dy- 
namical evolution the field x decreases. The function y(x) has one point of maximum y max = x m /2. 
If 2/max > then this evolution results in inflation (the region II) (if y max 3> 1> the evolution 
results instead in slow-rolling regime I). Then the field decreases and starting out from the y = v3 
we have the natural exit from the inflation (the region III) - inflation indeed ends spontaneously. 



Remark 9. The same types of reasonings can be conveyed for another cosmologically and quantum 
field theory popular model with the potential V = A</> 4 /4, where A > 0. Evidently, all the qualitative 
properties of its dynamics (inflation's spontaneous start and exit) will remain the same. Moreover, 
as follows from the analysis, the beginning of the inflation, its exit together with the presence of 
the slow-rolling regime doesn't depend on the model's physical parameters, be it either scalar par- 
ticle's mass m (for quadratic potential) or the coupling A for the 4 model. This demonstrate the 
universality of inflation, at least for the class of the power potentials. 
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Figure 3: The plot of numerical solution y(x) of the Eq. (fTTj) with initial data y(10) = 1.2. 



Remark 10. We have used above the analyticity of solutions of (|33|) at x > 0. This condi- 
tion will be of crucial importance in the studies on the (becoming very popular recently) effect of 
the crossing of the phantom divide line. In a series of works (cf., for example, [23], [25]) it has been 
shown that during the cosmic evolution it is possible to for kinetic term to become negative valued. 
In order for this to happen there should exist a region where kinetic term becomes equal to zero, 
i.e. (j) = 0. This implies the existence of a point x*, such that W(x*) = V(x*). Comparing this 
with ([ID]) it is possible to show that x* has to be a second order break point for the function y(x). 
Hence, for the discussed power potentials the effect of the crossing of the phantom divide line is out 
of the question. 

6. The Backlund auto-transformations for the Abel equation 

Throughout this paper we have used the Abel equation for construction of the general solution of 
the Einstein-Friedmann equations for the models of universe filled with scalar field with the given 
potentials. In this Section we'll show that one can use the established connection to study some 
interesting mathematical properties of the Abel equation itself. 

Theorem 3 The Abel equation Ul\) admits a set of Backlund auto-transformations: 

X {x)^x (l) {x), y(x)^y^(x), (42) 
satisfying the following properties: 

Property 1. If y(x,C) is a general solution of the ( T77]) with given x( x )> then applied to 
y and X generates a general solution 

x (x)^ X {1 \x), y(x,C)^y«(x,C), 

of a transformed equation 
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Property 2. This transformations can be applied arbitrarily many times; after the n-th step the 
n-times transformed function y( n > (x,C) will contain an arbitrary constant C, while x^ n ' = X ( x )- 

Property 3. The transformations can be inverted by setting the integration constant equal to 
zero. For example, for n = 2 

X(x) - X {2) = X {2) (x,CW=0) = X (x), 
where C^- 1 ' is a certain intermediate parameter. 

The proof of the Theorem 3. It was shown that using exact solution of the Abel equation (jlip 
y(x,C) it is possible to find the exact solutions x(t, to,C) (i.e. (f>(t,to,C)) and a(t,to,C) of the 
Einstein-Friedmann equations in the flat space-time. Define function ip(t,to,C) = a 3 (t,to,C). One 
can show that tp is the solution of the Schrodinger equation [8], [24] (with the notation ip = d 2 ip/dt 2 ): 

$ = 9tf (tty, (44) 

where U c (t) = V(x(t,C)). 

Remark 11. U c (t), unlike to V{x) does depend on constant C (for simplicity we'll omit con- 
stant to for everything that follows). In fact, using equation ([5]) it is possible to present t as a 
function of x and C: 

*-*-<»*khW£- (45> 

Inverting the (|45p (this can be done explicitly, as we will show below) we'll get x = x(t,C), i.e. 
V(x) = U c (t). This observation gives rise to what we'll call the reparametrization invariance 
of function U : 

U c (t) = U c (t(x,C)) = V(x). (46) 

This property will prove to be of utmost importance for the proof of the Property 1. 
Equation (JUJ) admits Darboux transformation [26], [27J: 



U a (t) - E/«(t) = U c (t) - I^Jog^C), (47) 



i,(t,c) -» ^(t,c,c^) = ^75y (i + / dt^ 2 (t,C) 



(48) 



where C^ 1 ^ is a new arbitrary constant of integration 0. Using the equations JT]) it is possible to 
define a new field variable 

,«=x«±2/«^log-l y , (49) 

as a function of time a^ 1 ) = x^^(t, C, C^ 1 )). Inverting this relation, we obtain t = t(x^,C^). Of 
course, the inversion of function in general is a highly nontrivial task, but the one easily accomplished 
in our particular case. 



7 Another one of new integration constants in (|48[) is introduced multiplicatively and cancels for all the formulas 
below 
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Proposition 5 The relationship x = x(t) can be explicitly inverted (as t=t(x)) for all explicitly 
known solutions of equation (or (G])j. 

In order to justify this statement one has to refer to (|45p . which is a result of inversion of x = x(t), 
defined by (fl9|) with the substitution tp^ — ► ip, Xq —> Xq. 

Using the relation t = t(x^\C,C^) it is possible to calculate the quantities U^}\t) = 
U^(t(xW,C,C^)). Assuming C« = C, one gets U£\t) = U^\t(x^ ,C)). Due to the 
reparametrization invariance (|46p : 

V^(x^) = U^(t(x^,C)). (50) 

I.e., the reparametrization invariance assures that a new potential V^' depends on the new inde- 
pendent variable x^ , but not on C . This is by all means true also for = 
On the next one is to calculate the quantity 

wW(xW,C) = (H^(t(x^,C)) X2 



where 

and ifj^ shall be defined from (|48[) with regards of the substitution = C . Finally, using the 
formula (fT0|) 

wW(xW,C) = V^(x^)e\y^(x^,C)), 
where 6 2 (y) is defined by (|31H one derives the thought-after " dressed" solution of Abel equation 



y (U (,(1) c) = , (51) 



Formulas (|50p and (|5ip defines the Backlund auto-transformations for the Abel equation (jlip and 
proves the Property 1. The proof of Property 2 is based on the usage of n-fold consequent Darboux 
transformations and the Crum determinant formulas [28j. The Property 3 can be checked directly: 
it is enough to let = (but before one lets this parameter be equal to C, because at this step 
C« and C are still being treated as the independent constants) in (|48j) and perform the second 
Darboux transformation: 

- tr® - l^iog^(t,c,cM = 0) = U a (t), 

Q.E.D. □ 
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